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Abstract 

We study conformal properties of the quantum kinetic equations in curved spacetime. 
A transformation law for the covariant Wigner function under conformal transformations 
of a spacetime is derived by using the formalism of tangent bundles. The conformal 
invariance of the quantum corrected Vlasov equation is proven. This provides a basis for 
generating new solutions of the quantum kinetic equations in the presence of gravitational 
and other external fields. We use our method to find explicit quantum corrections to 
the class of locally isotropic distributions, to which equilibrium distributions belong. 
We show that the quantum corrected stress-energy tensor for such distributions has, in 
general, a non-equilibrium structure. Local thermal equilibrium is possible in quantum 
systems only if an underlying spacetime is conformally static (not stationary). Possible 
applications of our results are discussed. 
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I Introduction 

In a regime where the curvature of a spacetime is less than the Planck scale, a semi- 
classical approximation to the quantum gravity is considered to be sufficient |]l|, |2|. This 
approximation is valid for "nearly" classical systems, for which quantum fluctuations in the 
stress-energy tensor are small enough In such cases the method of Wigner functions 
[Q, ^, ^ is particularly convenient. It allows one to derive, from quantum field equations, 
quantum kinetic equations for Wigner functions. Iterative solutions of the kinetic equations 
give quantum corrections to classical distribution functions, which are used to compute ex- 
pectation values of local observables such as the stress-energy tensor, the number-flux vector 
etc. The standard definition of the flat-space Wigner function exploits the Wigner- Weyl 
transformation |^] which fails in curved spacetime. Different approaches have been proposed 
to extend the notion of the Wigner- Weyl transformation to curved spacetime. In Ref. a 
Wigner- Weyl-type transformation is defined by considering the bundle of geodesic lines on 
a manifold. Ref. utilizes Riemann normal coordinate systems and adiabatic expansions 
similar to one used by |l2] for studying ultraviolet properties of Green functions in curved 



spacetime. In Ref. [10|, a covariant Wigner function is defined by using the formalism of 



tangent bundles. All three approaches mentioned, being intrinsically different, are locally 



equivalent to each other [11|, in the sense of adiabatic expansions, and yield the same quan- 
tum corrected Vlasov-type kinetic equations. The mathematical structure of the equations 
is rather complicated, even at lowest adiabatic order. This makes difficult the study of 
properties of generic solutions of the quantum kinetic equations in an arbitrary background 
spacetime. No doubt that any explicit computation of the Wigner function may be useful 
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for better understanding of quantum kinetic processes in gravitating systems. 

If a spacetime manifold is endued with some symmetries then both classical and quan- 
tum kinetic equations can often be solved explicitly. By now, only in a very few cases explicit 
solutions of the quantum corrected Vlasov equation in curved spacetime have been obtained. 



They include the Freidmann-Robertson-Walker cosmology [1^, 14], stationary spacetimes 



[15 1 and some others |16|. A common property of those solutions is that the quantum correc- 
tions are local in phase space. Once a solution of the Vlasov equation (a classical distribution 
function) is known, one can compute, to a given adiabatic order, the Wigner function which 
corresponds to this solution, by applying a certain finite order differential operator to the 
classical distribution function. This property implies that one can express the expectation 
value of the stress- energy tensor or of any other local observable in terms of moments of 
the classical distribution function multiplied by some local geometrical quantities. This fa- 
cilitates to a great extent the analysis of the back-reaction problem in those cases. 

The natural question arises how symmetries of a spacetime manifest themselves in 
Wigner functions. In the present paper, we consider a particular, but important, type of 
geometrical symmetries, the conformal symmetry. For simplicity, attention will be restricted 
to the case of a scalar field, spin-^ and spin-1 fields can be treated similarly [|l^, |l^]. The 
conformal symmetry is usually understood as the invariance of the generalized massless Klein- 
Gordon equation (with the term ^Rf involved Q) under the conformal rescaling of the metric 
and of the scalar field: gap — > a^Qap^f ^ja. If the scalar field is coupled to an external 
potential, the latter must also be transformed in such a way that the scalar field equation 
remains invariant under this extended group of conformal transformations. 

Two different views at the conformal symmetry are possible. When the external po- 



tential has a direct physical meaning, the conformal symmetry is a dynamical symmetry of 
the physical system under consideration. This is the case, for example, in the ip^ theory where 
the expectation value (ip^) plays a role of an external potential for quantum perturbations. 
The second point of view is formal. One can treat, for example, the mass of the scalar field 
as an external potential and use the conformal symmetry as a convenient tool for generating 



new solutions of the field equation |17]. 

In this paper, we study transformation properties of the covariant Wigner function 
under the conformal transformation. As we have mentioned, different curved-space exten- 
sions of the flat-space Wigner function are possible, all of them giving similar adiabatic 
expansions. The use of the tangent bundle technique is particularly convenient for exploring 
mathematical properties of the Wigner function. By using the formalism of tangent bundles 
we derive, to second adiabatic order, the transformation law for the Wigner function and 
prove the conformal invariance of the quantum-corrected Vlasov equation. We then use our 
result to find explicit quantum corrections to the class of distribution functions which are 
locally isotropic in momentum space, i.e. have the form Fci{x,p) = F(x, u°^{x)pa)- As it was 



shown by [18, 19 1, in the case where external fields other than gravitation are absent, such 
form of the distribution function imposes the following restriction on the underlying space- 
time: it must be or stationary or conformally static (with a special form of the conformal 
factor). In the presence of an external potential this condition has to be modified. We do 
this by using the conformal properties of the Vlasov equation. We also analyse the structure 
of the quantum corrected stress-energy tensor and number-flux vector, and show that the 
quantum corrections lead, in general, to non-vanishing heat flux and viscosity, even if the 
classical distribution implies local thermal equilibrium. 



The paper is organized as follows. In Sec. we outline the formalism of the Wigner 



function in curved spacetime. Section III derives the transformation law for the Wigner 



function. The lowest order quantum corrections to the isotropic distributions are evaluated 



in Sec. IV, and the quantum corrected stress-energy tensor and number-flux vector are 



calculated in Sec. |^ Section VI contains some remarks concerning possible applications of 



our results. In [Appendix A| , technical details used for the evaluation of the transformation 



law are presented, and in [Appendix B| , the conformal invariance of the quantum corrected 



Vlasov equation is proven. [Appendix C recalls some properties of the stress-energy tensor. 



Conventions used throughout are the following. Greek indices run from to 3, Latin 

indices from 1 to 3. The signature for the metric tensor is (H ). The Riemann tensor 

^^a(5 is defined by 

The Ricci tensor is Rap = R^avf5-> scalar curvature R = i?". 



[V„,V^]X'^ = < (1.1) 



II The covariant Wigner function 

To keep the paper self contained, we shall outline our definition and some properties 



of the covariant Wigner function in curved spacetime [10, 11|. Let <^(x) be a scalar field on a 
spacetime manifold M. with the metric ga/3{x), and Tx{M.) be the tangent space at a point 
X. The horizontal lift of the covariant derivative operator to the tangent bundle T{A4) 
is defined by [^] 

V„ = V. - ^ , (2.1) 
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where are components of a tangent vector pertaining to Tx{A4) , and T^^{x) are the 
ChristofFel symbols. Let us introduce the horizontal lift of the scalar field: 

*(x,y) =exp(y"V„)^(x) . (2.2) 

Then, the covariant Wigner function is defined as follows: 

fix,p) = {nh)-^ J^) f dSe-2^^>"/" {^ix,-y)^\x,y)) , (2.3) 

where the bracket stand for an ensemble averaging with some density matrix defined on a 
Cauchy hypersurface 0], the dagger indicates the hermitian conjugation. 

If the field f{x) is a solution of the generalized Klein-Gordon equation: 

h'^ V^V" -Ih'^ R{x) + V{x)] fix) = , (2.4) 



6 

then the Wigner function obeys the quantum corrected Liouville-Vlasov equation supple- 
mented by the quantum corrected mass-shell constraint 11]. The explicit form of the 



equations to second adiabatic order is given in Appendix B. We note here that a Wigner 
function satisfying these equations (to second adiabatic order) is represented as the series in 
derivatives of the Dirac J-function ||9|, PH : 

f{x,p) = [F,i{x,p) + n^Fgu{x,p)) 6{^)+n'Fi{x,p)6'{^)+n^F2{x,p)6"{^)+n^F^{x,p)5"'{^) , 

(2.5) 

with the argument of the 5-function being 

n = g''^{x)paP^ - V{x) . (2.6) 

The functions Fn{x,p) , n = 1,2,3, are expressed in terms of the classical distribution 
function Fci{x,p), in such a way that the quantum corrected mass-shell constraint is satisfied 



(see Appendix B), assuming that Fci{x,p) obeys the Vlasov equation: 



6{Q) [CF^{x,p)) = 0. 



(2.7) 



Here C is the LiouviUe- Vlasov operator |2^, 24 1: 



1 d 

2 dpo 



(2. 



Da being the horizontal lift of the covariant derivative operator to the cotangent bundle |2C]: 



d 

dpf3 



(2.9) 



The function Fqu{x,p) is found by integrating the quantum transport equation (see Appendix 



B) along classical trajectories 11 1 



Equation (2/7) defines the classical distribution function only on the mass shell 17 = 0. 
This is sufficient so far as one is interested in classical observables defined on the mass shell 



[22|. In quantum kinetic theory the situation is different |11]. Though the function Ff.i{x,p) 
is tied up with the (5-function in the expansion ( |2.5D , its properties out of the mass shell 
influence upon both the off-shell quantum corrections, the Fn's, and the on-shell quantum 



corrections, Fqu{x,p). The latter becomes evident if we get rid of the (5-function in Eq. (2.7) 
The function Fci{x,p) must satisfy the generalized Vlasov equation: 



£ Fci{x,p) = nAF{x,p) 



(2.10) 



Ap being an arbitrary function which is assumed to be non-singular on the mass shell. 
Then the equation governing the evolution of Fqu{x,p) explicitly involves the function A^? 



(see Appendix B| ). In the next section we shall see that, in general, one can not get rid of 
Ap. Even if the right-hand side of Eq. (2.10) happens to vanish for some system, for systems 



conformally related to the one under consideration, the Liouvihe-Vlasov operator annihilates 
Fci on the mass shell only. 

Ill The Wigner function and conformal transformations of a 
spacetime 

Let us now consider the conformal transformation on a spacetime manifold (we shall 
use throughout an over line for quantities in the conformally related spacetime): 

(p{x) = lp{x)/a{x) , (3-1) 

V{x) = V{x)/a{xf , 
a{x) being an arbitrary smooth function. Equation ( p.4| ) is known to be invariant under 
such a transformation Q. If (p{x) is a solution of Eq. (|2.4| ) on the manifold A4 with the 
metric Qapix) and in the presence of the external potential V{x), then Tp^x) is a solution of 
the conformally transformed equation on the manifold M with the metric gap{x) and in the 
presence of the external potential V{x), and vice versa. 

One can expect that the Wigner function f{x,p) associated with the field (p must 
be somehow connected to the Wigner function f{x,p) associated with the field ^. In the 
classical limit {h —>■ 0) the connection is very simple. The mass shell ( |2.6| ) is only scaled under 
the conformal transformation (we treat covariant components of the momentum vector, pa, 
and contravariant components of the vector as parameters which are not changed under 
conformal transformations, therefore, making difference between the tangent and cotangent 
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spaces). Thus, 

J7 = na~2. (3.2) 



Next, under the conformal transformation ( |3.1| ), the Liouville-Vlasov operator is transformed 
as follows: 



C = Ca + Cq, TT — O , (3.3) 
dpa 



where 



d 



By comparing Eqs. (p.2|),( p^ with Eqs. ( |B2D and (Bl), one can conclude that the transfor- 
mation law for the Wigner function should be 

f{x,p) = a(x)^ f{x,p) + quantum corrections . (3.5) 

We could also arrive at Eq. (^]^) by writing the transformation law for the Liouville- 
Vlasov operator in the form 

— d 

C = a~'^ C + naaT^ ■ (3.6) 
dpa 



From Eqs. ( ^.7] ) and (3^) it follows that the classical distribution function remains invariant 



under the conformal transformation. Then Eqs. (2.5) and ( |3.2| ) imply (3.5) 



One can also see from Eq. (p. 61) that, unlike Eq. (|2.7D , the off-shell Vlasov equation 
( 2.10| ) is not invariant under the conformal transformation. Namely, the function A^? is 



changed as follows: 



d 

Af{x,p) = Af{x,p) + aa-^ — Fci{x,p) . (3.7) 



In order to find the quantum corrections in Eq. ( |3.5| ) , we must turn to the definition of 
the covariant Wigner function (|2.3|). First, we write the transformation law for the operator 



Va , Eq. ( p.lD , when acting to a scalar function in the tangent bundle. Under the conformal 
transformation of the metric, Eq. ( |3.lD , it is transformed as follows (cf. Eq. ( B14D ): 

d 



(3.8) 



Let us now write the field ( p. 21) in the following way (compare Eq. (|3.l| )): 



V(a;) 



(3.9) 



The last factor on the right-hand-side of Eq. (3.9) is , by the definition, $(x, y). The second 
factor, which is an operator acting in the tangent bundle, can be expanded as follows (see 
Appendix A| ): 



1 



■Va + 



1 d 



3 ' " 2dy^ 



+ terms of higher orders , 



(3.10) 



where and i?" are given in Appendix A . 

Lastly, the first factor on the right-hand side of Eq. ( |3.9D is a function which involves 
the conformal factor a{x) and its derivatives. The known formula gives the expansion: 

A{x,y) := e^y-^-^a-'e^-y"^-^ 



a-' fl-y-a„ + i(y"a,)2-i 



a 



+ terms of higher orders , 



(3.11) 



where a^;/? := V/ja^. On the right-hand side of Eqs. ( 3.10| ), ( p. 11 ), we have kept only those 



terms which yield lowest order quantum corrections to Eq. (3.5). 



We are now prepared to write down the transformation law for the Wigner function, 



correct to second adiabatic order. Before doing that we note the following identity [11|: 
d 



y) = Vq, y) + terms of second or higher orders . (3.12) 



Then, with the aid of the results of Appendix A, we finally get: 



+ terms of higher adiabatic order , (3.13) 



where the operator Da is conformally related to ( |2.9D , and A"' (dp) and B°^{dp) are the 
following operators (cf. Eqs. (A2),(A3)): 

A-{dp) = (2a^5^ - a"g^J , (3.14) 

B^{dp) = {2a^,, 6^ - a'^g,^ - 8a^ a, 5" + 2a^ a^^, 6^ + 4a- a^g,^) Q^Jp^g^^ ■ (3-15) 

We have denoted Uf^-i, := Vj^a^ and a" := 'g°^^ap. 

Equations (|3l^)-(|3l^) give the transformation law for the Wigner function which 



we have been seeking, correct to second adiabatic order. It is worth noting that the transfor- 
mation law has been derived irrespective of any specific form of quantum kinetic equations, 
and it holds independently of whether or not the Wigner functions satisfy these equations. 
One could follow a different way. Given equations which govern the evolution of a Wigner 
function, one may seek a transformation which leaves the equations invariant under the con- 
formal transformation of the metric and the external potential. This is what we did in order 
to obtain the transformation law at zeroth order, Eq. ( |3.5| ). However, this program meets 
serious problems at higher orders because the structure of the quantum kinetic equations in 
curved spacetime is very complicated. Instead in [Appendix B we prove, to second order, the 
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conformal invariance of the collisionless quantum kinetic equations by using the transforma- 
tion law for the Wigner function. 

One may make sure that Eqs. (|3.13| )"( |3.15| ) yield the right transformation law for 



physical observables. For example, the number-flux vector defined by |11| 



—==Paf{x,p), (3.16) 
^-9\.x) 

is transformed as follows: 

(-5)^/^Ja) = (-5)'/'(Ja). (3.17) 

The transformation law ( ^.17 ) implies that the number of "particles" (defined as in kinetic 



theory) in a comoving three-volume, dN := {J°')dT,a, is an invariant of the group of conformal 
transformations. Here dSo, := Eacr^iv dx" dx^dx^ . 

The transformation law for the stress-energy tensor depends upon an explicit form 
of the external potential. For a scalar field coupled to curvature one gets (see Appendix Cj ) 



{g/gf/^ (T„^) = (T„^) + ^ (6^ - 1) (2aa(, V^) - ag^^a^V" -g^pUl) {N/a) , (3.18) 



where ^ is the nonminimal gravitational coupling constant and 



N{x)= I ^Pf{x,p). (3.19) 



IV Quantum corrections to isotropic distributions 

Let us consider a system in external fields ga/3 and V whose classical distribution 
function (which is assumed obeying the Vlasov equation) is locally isotropic in momentum 
space, that is Fci{x,p) = F{x, u°'{x)pa), u"(x) being some world velocity field. Of course, the 
external fields must be consistent with the symmetry of the distribution function. From the 
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results of Refs. |jT^, |T9| and from the conformal invariance of the Vlasov equation it follows 
that only two possibilities exist: 

A. The metric is conformally stationary, i. e. the line element takes the form (in preferred 
coordinates) 

ds^ = a^(dt^ + 2%idx°dx' + g^jdz'-dx^) , (4.1) 

with '9ij being functions of the spatial coordinates alone, and either the potential 

V equals zero or the function V := Va? does not depend on time; 

B. The metric is conformally static, that is = in some coordinate system, and the 
function V depends on time alone (since the Vlasov equation is linear in V , one can conclude 



that in Case B the function V may be of the form |19| V = Vi{t) + V2{x^) ] we shall not 
deal with this more general situation here). 

In accordance with our strategy, we can solve the quantum kinetic equation in an 
auxiliary spacetime with the line element 

df = + 2gQidx°dx' + g^^dx^dx^ , (4.2) 

and then use Eq. ( 3.13| ) to calculate the Wigner function in the physical spacetime. 

— a 

In the spacetime with the metric 'g^p there exist a unit time-like Killing vector ^ , 



^ = 1 (^ =5^111 the preferred coordinates) , which obeys the equation |g3|] : 

Va?/3 + Vpla = in Case A , (4.3) 

or Vc^p = in Case B . (4.4) 

The potential must satisfy the conditions: 

TV;a = ^ = Ttv,^p = ... in Case A, (4.5) 

A"^F;„ = = A°^F;„/3 = • • • in Case B . (4.6) 
12 



Here 

A°^:=rf-r^ (4.7) 

is the space-like projection operator and V-a ■= etc. 

In the two cases considered, the locally isotropic distribution functions take the 

form (one can see that in Case B the distribution function bellow has the desired form 
F(x,u'^{x)pa) on the mass shell): 

Fci=F{Tpa) in Case A, (4.8) 



F,, = F(Va"Vp/3) in Case B. (4.9) 



We proceed now to the derivation of the quantum corrections to the classical distri- 



bution functions. Both the function ( [4. j ) and the function (O) satisfy the off-shell Vlasov 
equation in the external fields 'g^p and V, with Ap being zero (see Eq. (^loD). Note, how- 
ever, that for systems conformally related to the ones under consideration, Ap ^ 0, as it 
follows from Eq. (p.7[). 

The semiclassical Wigner functions which correspond to the distribution functions 
(|4.8| ) and (^^) in the spacetime M are represented by the same expansions as in Eq. ( |2.5[) , 
with the argument of the (5-function being 

^ = r^PaPp - V . (4.10) 

The off-shell quantum corrections are obtained by substituting the classical distribution 
functions ( [4.8| ) and ( |4.9| ) into Eq. (|B6|). Both in Case A and in Case B, they can be written 
in the form: 
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T3 = ^{2V^P„P^-V"V„) F^,. (4.11) 
Next, the on-shell quantum corrections are found from Eq. ( |B7] ). It follows from the 



conditions listed in Eqs. (4^) and ( [4 .61) that the equation for Fqu can be integrated as if 



the potential V were constant. A solution of this equation which corresponds to the classical 



distribution ( [4.S| ) was found in Ref. [15|. It reads (with C ?a = 1 



_ 1- 92 

Fqu = —^ Ra/3 g^-g^ ^ci m Case A. (4.12) 



Let us now consider Case B. It follows from Eq. ( [4.4] ) that the following identities 
hold |2|: 

C Ra/Sfiu = = ^ Ral3fiu;a ■ (4-13) 

By using Eqs. ([4. 4]) , (|4.6|) and ( [4.13 ) one readily gets that the function 



_ 1 _ 92 

F'lu = -t:^ Rap Fa in Case B (4.14) 

12 dpadpp 



satisfies the quantum kinetic equation ( |B7| ) , with being of the form 

Equations (|4.11| ) and ( 4.12 ) give the second order quantum corrections to the locally 



isotropic distribution (|4.8| ) in a stationary spacetime in the presence of a static potential 
(Case A), while Eqs. ( [4.11| ) and ( [4.14| ) give the quantum corrections to the distribution ( [4. 9] ) 
in a static spacetime in the presence of a homogeneous potential (Case B). As we have said, 
the Wigner functions for systems conformally related to the ones which have been considered 
in this section are obtained by using Eqs. ( |3.1| ) and ( |3.13| ). In the next section we shall 
analyse the structure of the quantum corrected number-flux vector and stress-energy tensor 
in these cases. 
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V Quantum corrections to a perfect fluid 



The distributions ( [4.^ ) and ( |4.9| ) both imply that the classical stress-energy tensor 
has a perfect fluid structure [p!9| ], 

(T^f ) = (e,i + P,i) r f - , (5.1) 

the classical energy density Ed and pressure Pd being given by 

Ed = 27T / dpp^ (V+p^)'/^ Fd{V,p) , 
Jo 

Pd = -2Tr dpp^{V+pY^/^Fd{V,p), (5.2) 
o Jo 

with (compare Eqs. (gj), (|49|)) 



_ F{JV + p2) in Case A 

Fd{V,p) = { (5.3) 

F{p) in Case B . 

The classical number-flux vector associated with the locally isotropic distributions reads: 

(Jd>="c/r, (5.4) 

the classical number density being 

/•OO 

rid = 27, dpp^FdiV,p). (5.5) 

JO 

It should be mentioned that we have taken into account positive "energies" only, 



E = \JV + p^, when integrating over the mass shell. This is tantamount to multiplying the 
distribution function by the step function ||2^ 0(C"Pa). Though ^"'pa is not a constant of 
motion in Case B, the step function is "ignored" by the Liouville-Vlasov operator, neither 
it contributes to the quantum kinetic equation, since differentiation of the step function 
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gives which vanishes on the mass sheh (we assume that, in the classical limit, V is 



non-negative everywhere in the spacetime). Note that negative "energies", E = —yV + p"^, 
correspond to antiparticles |0] . 



The quantum corrections to Eqs. (|5.lD , (|5.4D are easily obtained by substituting the 



Wigner functions found in the preceding section into Eqs. (C4) and (|3.16| ). We shall first 



consider Case B. The second order quantum corrections to the number flux ( p. 41) can be 
shown to vanish. The quantum corrected stress-energy tensor can be written, after some 
simple algebra, in the following form: 



^^—2_a0 /3 , 1, 1 1 _ \ 1 —2-a-B 

-h'V A (^-(e--)Mi,3--M2,4--V^Mi,4j+-F e rMi,3. ( 
Here V := etc. , A"^ is defined by Eq. (O), and the M's denote the integrals: 



Jo 



The stress-energy tensor ( ^.Gj) has an "almost" perfect fluid form, except the second 
term which involves R . The latter has in general off-diagonal components. Note, however, 
that the Einstein equations impose severe restrictions on the underlying spacetime. Suppose, 



for example, that the conformal factor in Eq. (|4.l| ) depends on time only: a = a{t) (recall 
that we are considering Case B, that is ggi = in the preferred coordinates). If the total 
stress-energy tensor of classical matter is assumed to have a perfect fluid form then the 



physical spacetime is necessary a Robertson-Walker one 11. Then R = 2K A , with 



K = —1,0, 1, that is the quantum-corrected stress-energy tensor preserves the perfect fluid 
structure in this case. 

Let us write down the stress-energy tensor with the second order quantum corrections 
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in a Robertson- Walker spacetime with the hne element 

ds^ = a' {dt' - - r^d9^ - sin^ 9d^^) . (5.8) 

1 — A 

Since the Bianchi identities hold automatically for the stress-energy tensor derived in kinetic 
theory, it suffices to know the energy density e := (TqO). Making use of Eqs. (U) and 

(|5.6|) gives 



+ ^(1 - Q^)^{U + 2-U)M,,2 + ^{U + 2-UfM,,s. (5.9) 

Here dots stand for the time derivative, and the M's are given by Eq. (|5.7l ), with V = Ua?, 
U being the classical part of the external potential in the Robertson- Walker spacetime (see 
Appendix C| ). Equation ( |5.9| ) extends the result of Ref. [13| to the case of a variable potential 



and the nonminimal coupling to curvature. Physical manifestations of the quantum kinetic 



corrections in the Friedmann-Robertson- Walker cosmology have been examined in Ref. |14]. 
It was shown there that the quantum corrected energy density causes an apparent time 
variation of the value of the gravitational constant in the early universe. 

We turn now to Case A of the preceding section. The structure of the quantum 
corrections is more complicated in this case. To avoid lengthy expressions, we shall confine 
our attention to systems for which V = 0. Substituting the distribution function (|4.8| ) and 
the quantum corrections to it given by Eqs. ( [4.1l| ), (|4.12| ) into Eqs. ( |3.16D and (C4) gives the 



quantum corrected number-flux vector and stress-energy tensor in a stationary spacetime 
(recall that ^"^^ = 1): 
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Here J^^^^ is the Riemann tensor of the stationary spacetime, and 



(5.11) 



Mn = 2TT / dpp^ F{p) . 

Jo 



(5.12) 



The remarkable feature of the quantum corrections in Eqs. ( p. 10 ), ( 5.11| ) is that they 
have a non-equihbrium structure. To show this exphcitly, let us compute the invariants of 
the number-flux vector and stress-energy tensor. The hydro dynamical velocity associated 
with the number flux (15. 101) is 



(5.13) 



If we decompose the stress-energy tensor ( p. 11 ) like in Eq. ( |C8| ) we find next, with the aid 
of Eq. ( piC| ), that the eigenvector of (T"^) is 



(5.14) 



In general. If does not coincide with v". This implies that the heat flux defined by ||22| , 



(e + P) (u" - U°) 



(5.15) 



does not vanish. Indeed, substituting Eqs. ( p.l3D ,( 5T4| ) into ( 5.1E| ) gives 



r = — (2M0M3/M2 - Ml) R^u 



(5.16) 



The combination of the moments on the right-hand side of Eq. (5.16) can be zero for a very 
special choice of the function F{p). It is surprising that the heat flux as well as the viscosity 
(see Appendix C| ) do not vanish even when the classical distribution function describes local 
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thermal equilibrium ||2^. Though the effect is purely quantum, it could play an important 
role in the early universe and in the evolution of massive stars. Some physical consequences 
of this will be considered elsewhere. 

If (and only if) the manifold Ai is static, that is the classical fluid motion is irrota- 
tional, the quantum corrections in Eqs. (|5.10| ),( p.llD vanish, which is consistent with Case 
B in the limit V ^ 0, ^ = 1/6 (compare Eq. (^.6|). 

Finally, the quantum corrected number-flux vector and stress-energy tensor in a con- 
formally stationary spacetime can be easily found with the aid of Eqs. ( 3.17] ),( p.l8| ). For 



the conformal coupling (^ = 1/6) they are obtained by simply multiplying Eqs. ( |5.10| ) and 
( 5.11| ), respectively, by a""^ and a~^. 



VI Concluding remarks 

In summary, we have considered conformal properties of the covariant Wigner func- 
tion. The transformation law ( |3.13| ) relates Wigner functions in manifolds conformally related 
to each other, correct to second adiabatic order. Given a solution of quantum kinetic equa- 
tions which govern the evolution of a quantum distribution function (a Wigner function) 
of a specific system, one is able to evaluate Wigner functions for a wide class of systems 
conformally related to the one under consideration. 

There exist two possible applications of the result. First, our method allows one to 
simplify an original system by reducing it to a simpler one for which the evolution can be 
easily solved. As an example, we have found explicit solutions of the collisionless quantum 
kinetic equations in conformally stationary /static spacetimes and analyzed the structure of 
the quantum corrected stress-energy tensor and number-flux vector. 
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Second, in some theories of gravity, like Brans-Dicke-Jordan-type ones which have 
been attracting a great deal of attention in recent years in connection with the low energy 
limit of string theory (see, for instance, Ref. [^), two metrics are involved: a physical metric 
which enters matter equations, and the so called Einstein metric which satisfies Einstein- 
like equations. The two metrics are usually conformally related to each other (though more 
general relations have been suggested, see Ref. [^]), and one might want to solve matter 
equations on the physical manifold and then express the solutions in terms of the metric 
associated with gravity. Our method described in the present paper is especially suitable for 
such theories. 

We would like to finish our discussion by the following remark. It is well known 
18, 19, 22, 24] that local thermal equilibrium in a relativistic gas is possible in a very restric- 



tive class of physical spacetimes. As we have shown, in a quantum system thermal equilibrium 
is even more exceptional than in a classical one. Spacetimes in which quantum fields can 
be in local thermal equilibrium are restricted to conformally static ones. Although we have 



analyzed collisionless equations only, this is also correct for collision-dominated systems |24], 
because quantum corrections to the collision integral are supposed to be small compare with 
ones to the Liouville-Vlasov operator when a system is nearly equilibrium. Thus our results 
provide further motivation for studying non-equilibrium processes in quantum systems in 
curved spacetime. 
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Appendix A : Evaluation of the transformation law for the 
Wigner function 



One can derive Eq. ( 3.13| ) by applying the Campbell-Baker-Hausdorff formula |2S]: 



expiD) exp{-D)=exp(^D-D-^[D,D] + ^[D + D,[D,D]] + ..)j , (Al) 

for the operators y'^Va and y'^V a- Making use of Eq. (|3.8D,( |AT| ) and taking into account 
the fact [0] that the operator Vq, commutes with and annihilates immediately give 
Eq. (|3lc| ) with 

= 2y'^y-a, - y^y.a^ , (A2) 
= y'^VaA'' - 

a.. (A3) 
= 2y-y^y'^(V^a,) - y^^y'yAV^a^) - 8y"{y''a,f + 2y^y-'y,a^'a^ + 4.y-'y,y^a^a" . 
Next, with the aid of Eqs. (|3.10D and we obtain: 



Z{x,y)^{x,y) = (^l-i^^V^ + ^S^^ + ^A^A'^V^V^) *(x,y) 

+ terms of third or higher orders . (^4) 



Using Eqs. (|3^,(|3ll|),(|A^) and analo gous identities for the conjugate field, we arrive at 



-y)^Hx,y) = a-^l- y'^y^a^.p - ^-A^V^ - ^ 



2 6 

+ ^74"^^VqV^ j —y)^\x, y)j + terms of higher orders . (A5) 



Finally, substituting Eq. ( |A5| ) into the definition of the Wigner function ( |2.3D and integrating 
by parts give Eq. ( p.l3D (the last term in the bracket on the right-hand side of Eq. (A5) 
leads to terms of fourth adiabatic order which we do not keep). 
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Appendix B : Conformal invariance of the collisionless quan- 



tum kinetic equations 



The equations for a Wigner function of a scalar field obeying the generalized Klein- 
Gordon equation ( |2.4| ) read ||8|, |l^] , to second adiabatic order: 



cf{x,p) = n'kf{x,p) , 
nf{x,p) = -h^uf{x,p) , 



(Bl) 
(B2) 



where Q and C are respectively given by Eqs. (2.6) and (|2.^), and the operators on the 



right-hand side of Eqs. (|Blj) and B2 



are 



Q3 



1 92 1 

A = -RufSf.aP'" - —Rc,f„3u-^„p^'p''- 

6 '^'^ OPaOp/s 24 ^'^ dpaOppdpo 



12 "dpa 



1 ^2 
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ril:ryT^ h — V-al3a-- 



Q3 



dpadpp 24 ' dpa 48 ' dpadppdpa 



1„ ,.5 1.. 92 



' 8' '"^'dpadpfs 



(B3) 



(B4) 



Here Da is defined by Eq. (12. 9|) , semicolons signify the covariant differentiation associated 
with the metric Qfj^iy. 

The operators A and 11 satisfy the following identity 



(B5) 



which allows one to look for solutions to Eqs. (^, (|B2|) of the form |, The F's 
in Eq. (p. 51) are given by 



Fi(x,p) = IlFci{x,p) 
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Fsix,p) = ^[[U,Q],n]Fa{x,p) , (B6) 
while the function satisfies the equation [pl|]: 

C Fgu{x,p) = A Fa{x,p) + 11 Af{x,p) , (B7) 



Ap{x,p) being defined by Eq. ( 2.10| ) 



Consider now a conformally related spacetime with the metric g^f^ = Qap /o? and sup- 
pose that f{x,p) satisfies the quantum corrected Vlasov equation and mass-shell constraint 



in the presence of the external fields g^R and V = Va^. As we have shown in Sec. ID, 



/a/3 

the conformal transformation of the spacetime induces the following transformation of the 
Wigner function defined in the cotangent bundle: 

fix,p)=a^il + n^C)7ix,p) , (B8) 



where the operator C is, to second adiabatic order (cf. Eq. (|3.13D , 



C = — ( -4a^;v + a°'^g^^ + 28a^a„ - 10a" Uag^^ + ^a^D^ - 3g^^a"Dc 



l/_ __ __ __\<9^ 

+ ^ [-2a^.yPa + a°'„g^^pa + ^a^a^Pa - 2a°'aag^tyPa - Ag^^aaa^p^ j ■ (^^) 

If one takes into account the transformation laws for the Wigner function, Eq. ( |B8D , 
and for the Liouville- Vlasov operator, Eq. ( |3.3| ), and uses the equations for the function 
f{x,p), one arrives at the following equations: 

^f{x,p) = h^(t+[±,c]-ao,^ifi+[c,n]))7ix,p) 

+ terms of higher orders , (BIO) 
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nfix,p) = -h'{u + [c,n]] fix,p) 



+ terms of higher orders . 



(Bll) 



Now, our task is to prove that these equations are tantamount, to second adiabatic order, to 



Eqs. (Bl), (B2). Tedious but straightforward calculations give the following results: 



— - — - d — 1 9^ 



dpo 



(B12) 



A+[C,C] = (A + a„^n)a2 + i(a^., + a°a„g^J^^^£ 



i_ _ __ __ - 

- — (a;.;i.a-8a^;^a<, + 4a°a^;a-6a°o«a^(7^^)^^P^^-^Sl. (B13) 



When deriving Eqs. (B12), ( B13| ) we have used the transformation laws for the Christoffel 
symbols and Riemann tensor [23|: 



^ UI/ ~ ^ 111/ ~ ^ LL^U ~^ ~ QnU^ ? 



(B14) 



O '^RafiPu — Ra/iPu — '^0'v[a9^](3 — '^(^ I3[a9 



(B15) 



where 0^/3 := Oq,;^ — 0^0^ + \ 'g^pa^ay. Recall also that Oq, := 9a In a, a" := 'g^^ap and 
semicolons in Eqs. ( [B9D , ( |B12 )-( B15 ) signify covariant differentiation associated with the 
metric Qap- 

Substituting Eqs. (^), (^) into Eqs. (^), j W^ ) completes the proof. 



Appendix C : The stress— energy tensor 



Consider, for simplicity, the theory of a real scalar field described by the action: 



S = ^ I (fx 



(CI) 
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R being the Ricci scalar, U an external potential not involving a metric dependence. 

The variation of the action with respect to the metric tensor yields the stress- energy 

tensor 



We have introduced, for convenience, the generalized external potential |28] 

V = U-n\C-]:)R. (C3) 
o 

The expectation value of the stress-energy tensor can be expressed in terms of the covariant 



Wigner function (|2.3|) as follows |11]: 



(T^,(e, V)) = f /'^ . (p^p. + h'^R^u) f{x,p) 

J \/ — Q{X] 



V-9{x) 

+ ^'(7 - 0(V;.V, - / f{x,p) . (C4) 

4 J ^/-g{x) 



The action (|C1| ) is invariant under the conformal transformation (3J.). Notice that 
the potential U is transformed as the following: 

U = U/a^ + Ti^ (1 - 60 a"^ V^a . (C5) 
Since the transformation law for the potential explicitly involves the dependence on the 



metric, the stress-energy tensor is not conformally invariant. With the aid of Eq. ( 3.13 ) we 
obtain the following relation: 

(6^-1) (^a(^V^) - ^ g^^ a„ V° - a^a^ - ^ g^^ (V"a«)^ j f-^. 7(^' P) ■ (C6) 
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Simple algebra then leads to Eq. ( p. 18 ). When deriving Eq. (C6) we have used the identity 



(B15) and the following property of the operator ( |2.9D : 

f^D^ f{x,p) = V, / ^ fix,p) , (C7) 



which was proven in Ref. |11]. 



From Eq. ( |C(]| ) it follows that the stress-energy tensor (multiplied by the density 
(—5)^/^) is conformally invariant, but not tracefree, if ^ = 1/6. vanishes only if y = 0. 
Any physical stress-energy tensor can be decomposed as follows 

Tf,^ = {e + P)u^u^ - Pg^^ + U^^ , (C8) 

where £ and P are, respectively, the energy density and the isotropic pressure, is the 
tracefree viscosity tensor [H" = = u^'Hap], and Ui, is a unit time-like vector [u'^Ua = 1]- 
The perfect fluid stress-energy tensor takes the form: 

,0 

T^,u= (e + P) u^u^ - P Qf^u ■ (C9) 

1 

For small deviations from the perfect fluid structure, T^i, =T^u + Tf_iu, one easily finds, to 



first order in the deviations: 



o 1 

e +u°'u'^ Tap 

P +iA"'3 Ta^ 



(CIO) 



u 



Tap /{e + P) 



where the projection operator has been introduced: 



^i^v =u^u^ -g^u ■ (Cll) 
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